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SUMMARY 


An  optimization  design  tool  has  been  developed  for  a  flywheel  rotor  of  fiber- 
reinforced  composite  materials.  This  tool  specially  considered  the  process-induced 
residual  strains  using  an  effective  temperature.  In  the  second  period  of  the  current  project, 
most  endeavors  have  been  devoted  to  measure  the  residual  strains  of  both  helically 
wound  rotors  and  hybrid  multi-ring  rotors.  We  have  developed  a  new  epoxy  system 
which  can  decrease  the  residual  strains  by  about  60  %  compared  to  the  conventional  high 
temperature  cure  epoxy  system.  This  report  consists  of  three  parts:  Part  I.  The 
investigation  of  new  epoxy  systems;  Part  II.  The  effects  of  the  helical  winding  angles  on 
the  residual  strains;  Part  III.  Development  of  flywheel  rotor  optimization  tool,  FlyOpt 
considering  the  residual  strains. 

The  first  part  discusses  the  effects  of  filament  winding  angles  on  the  residual  strains 
and  through-thickness  material  properties  of  a  thick  filament  wound  composite  cylinders. 
The  filament  wound  cylinders  of  37.5  mm  inner  radius  with  a  75  mm  radial  thickness 
were  manufactured  using  E-Glass  fibers  and  Epoxy  resin.  The  cylinders  were  sliced  along 
the  radial  direction,  and  the  through-thickness  properties  and  the  residual  strains  were 
measured.  The  through  thickness  properties  were  measured  from  90  (hoop)  to  60  degree 
wound  cylinders,  i.e.,  the  radial  Young’s  Moduli,  the  coefficients  of  thermal  expansion, 
and  the  radial  tensile  strength.  Significant  decrease  in  the  residual  strains  in  the  radial 
direction  was  noticed  as  the  winding  angles  decrease;  1.23,  0.99,  and  0.54  X  10'3  for  90, 
75  and  60  degree  winding  angles.  The  residual  strains  were  successfully  predicted  using 
the  effective  temperature  drop  of  140  °C  in  both  two-dimensional  analysis  and  three- 
dimensional  FEM. 

The  second  part  discusses  the  effects  of  rotor  sizes  and  different  epoxy  systems  on 
the  process-induced  residual  strains  accumulated  within  wet-filament  wound  thick  hybrid 
multi-ring  composite  rotors.  Three  different  epoxy  systems  which  cure  at  different 
temperature  profiles  were  considered:  high,  middle,  low  (room)  temperature  cure  epoxy 
systems.  To  evaluate  the  characteristics  of  the  epoxy  systems,  the  radial  stiffness  and 
strength  and  the  coefficients  of  thermal  expansion  were  measured  from  the  Glass  and 
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Graphite  rotors.  Hybrid  two-ring  composite  rotors  of  Glass  and  Graphite  fibers  with 
equal  thickness  were  manufactured  for  three  different  inner  to  outer  radius  ratios,  i.e.,  1.5, 
2.0  and  2.5.  A  split  ring  method  was  then  used  to  measure  the  radial  and  circumferential 
residual  strains  within  the  rotors,  following  liquid  penetration  tests  to  detect 
delaminations  on  the  surfaces.  Delaminations  were  observed  in  the  high  temperature 
cured  rotors  of  radius  ratios  2.0  and  2.5.  The  maximum  radial  stresses  are  observed  to 
occur  always  at  the  interface  between  the  Glass  and  Graphite  rings.  The  measured 
residual  strains  are  successfully  predicted  using  a  two  dimensional  stress  analysis  with  an 
effective  temperature  drop.  It  is  shown  that  the  middle  temperature  cured  rotor  attains  as 
high  strength  as  the  high  temperature  cured  rotor  while  reducing  significantly  the  residual 
strains. 

The  third  or  final  part  present  how  the  stress  analysis  module  was  developed  and 
linked  to  an  optimization  program  in  which  the  design  variables,  e.g.,  the  thickness  of 
each  hybrid  ring  and  winding  angles  can  be  optimized  for  the  various  size  and  material 
sequences  of  the  composite  flywheel  rotor  with  a  consideration  of  both  the  residual 
stresses  and  the  rotational  centrifugal  stresses.  The  program  FlyOpt  has  been  further 
updated  to  incorporate  the  current  experimental  results  and  to  have  a  user-interface 
convenient  for  designers  of  a  flywheel  composite  rotor. 

Currently,  the  developed  program  is  successfully  utilized  by  Structures  and 
Composites  Laboratory  in  Stanford  University  headed  by  Prof.  Steve  Tsai,  and  Korean 
Institutes  of  Machine  and  Machinery  in  Daedug,  Korea.  As  a  result  of  the  current 
research,  two  papers  were  submitted  for  technical  publication:  the  manuscript  in  the  Part  I 
to  Composites  Science  and  Technology,  and  the  other  manuscript  in  the  Part  II  to  Journal 
of  Composite  Materials. 
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PART  I.  EFFECTS  OF  WINDING  ANGLES  ON  RESIDUAL  STRAINS 
AND  THROUGH-THICKNESS  PROPERTIES  OF  THICK 
FILAMENT  WOUND  CYLINDERS 

1.  INTRODUCTION 

Filament  winding  has  matured  to  become  one  of  major  manufacturing  process  in 
the  fabrication  of  high  performance  fiber  reinforced  composites.  As  applications  of 
polymer  composites  increase,  there  has  been  a  great  demand  for  thick  sections  of  fiber 
reinforced  composites  such  as  flywheel  rotors  for  energy  storage  system,  filament  wound 
hulls  for  submersibles,  and  monocoque  hulls  for  armoured  fighting  vehicles  [1—3]. 
However,  such  thick  sectioned  structures  inevitably  generate  residual  stresses  which 
mostly  deteriorate  the  structural  performance  [4,5].  Attempts  should  be  thus  made  not 
only  to  decrease  such  residual  strains  but  also  to  predict  them  in  a  design  process  [6]. 

Most  design  for  filament  wound  structures  was  based  on  netting  analysis  which 
uses  the  in-plane  material  properties.  It  was  simple  and  widely  accepted  particularly  for 
relatively  thin  sectioned  structures.  However,  as  the  thickness  increase,  a  more 
sophisticated  analysis  such  as  finite  element  methods  is  needed,  and  through-thickness 
material  properties  should  be  defined.  Recently,  attempts  were  made  to  determine  the 
through-thickness  material  properties  of  FRP  materials  [7,8].  In  this  work,  compressive 
and  tensile  tests  were  performed  on  various  perform  types  such  as  unidirectional  prepreg, 
continuous  hoop  wound  filaments,  0/90  woven  prepreg.  However,  helically  wound 
structures  were  not  tested  for  various  winding  angles.  Coefficients  of  thermal  expansion 
(CTE)  of  the  helical  wound  structure  in  the  through-thickness  direction  were  never 
measured:  this  material  property  is  also  very  important  particularly  for  analysis  of 
residual  strains  within  the  thick  wound  structures. 

Prediction  of  the  residual  stresses  generated  during  the  curing  process  in  thick 
composite  cylinders  generally  requires  an  understanding  of  very  complicated  multi¬ 
physics  including  a  detailed  universal  rheological  model,  the  process  of  gelation,  cure 
kinetics,  the  exotherm  and  interaction  between  the  cylinder  and  the  mandrel.  Several 
investigators  have  developed  theoretical  models  of  the  filament  winding  and  curing 
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processes  [4,5,9-12].  In  attempts  to  reduce  the  residual  strains  during  filament  winding 
process,  a  multi-stage  cure,  or  an  in-situ  cure  and  a  continuous  curing  process  can  realize 
the  residual  stress  advantages  over  the  conventional  one-stage  cure  [13,14],  Nevertheless, 
the  one-stage  cure  has  been  still  widely  used  in  the  composite  material  industry,  and  it 
can  be  more  productive  as  long  as  the  residual  stresses  are  properly  managed.  The 
residual  strains  are  mostly  caused  by  the  temperature  drop  during  the  cool  down  process 
and  the  chemical  shrinkage  experienced  within  the  matrix.  An  effective  temperature  drop 
which  amounts  to  the  residual  strains  was  determined  by  measuring  the  residual  strains 
and  using  two-dimensional  stress  analysis  [6]. 

In  this  study,  the  effects  of  the  helical  winding  angles  on  the  through-thickness 
material  properties  and  the  residual  strains  of  the  helically  wound  structures  are  studied. 
Two-dimensional  stress  analysis  with  an  effective  temperature  drop  will  be  used  to 
estimate  the  stresses  and  strains  accumulated  inside  the  helically  wound  structures.  The 
through-thickness  material  properties  such  as  radial  stiffness  and  strength,  and  the 
coefficients  of  thermal  expansion  (CTE)  for  different  helical  winding  angles  will  be 
measured  and  used  in  the  analysis  [15].  A  destructive  split-ring  test  is  then  performed  to 
measure  the  residual  strains  within  the  helically  wound  structures.  The  measured  residual 
strains  will  be  compared  with  results  from  the  two-dimensional  analysis  and  three 
dimensional  finite  element  methods. 


2.  STRESS  AND  STRENGTH  ANALYSIS 


The  adjacent  layups  having  the  alternating  angles  [±<j>\  act  as  orthotropic  units,  and 
the  helically  wound  cylinder  is  assumed  to  be  cylindrically  ortho  tropic.  A  constant 
orientation  angle  <f>  throughout  the  cylinder  is  considered,  and  the  stress-strain- 
temperature  relations  for  an  axisymmetric  multi-layered  cylinder  in  cylindrical  coordinate 
are  given  as  Figure  1 . 
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where  ar,  cr&  sr,  eg  are  the  radial  and  circumferential  stresses  and  strains,  respectively. 
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The  components  of  the  compliance  matrix  Sy  are  obtained  by  transforming  material 
properties  from  the  material  axis  to  the  cylindrical  coordinates  [16,17].  The  axial  stresses 
are  assumed  to  be  zero  based  upon  the  plane  stress  assumption.  The  coefficients  of 
thermal  expansion  (CTE)  ar  and  ccq  are  also  obtained  from  material  properties  defined  in 
the  material  axis.  An  effective  temperature  change  ATe  will  be  determined  to  account  for 
the  residual  strains  generated  in  the  cured  rotors.  The  residual  strains  are  mostly  caused 
by  the  cool  down  process  and  the  matrix  chemical  shrinkage. 

The  equilibrium  and  compatibility  equation  together  with  the  stress-strain 
relationship  in  Eq.(l)  yield  the  radial  and  circumferential  strains  and  stresses  for  a 
cylinder  [16-18]: 


£r  =  kC^7k~'  +  kC2(ps~k~x  +  <Pt2A^ 

(2a) 

£g  =  C]q>]7K~'  -  C2(p-/~K~X  +<pT2 ATe 

(2b) 

a0  =  kC/‘ -  kC27^]  +  q>nA  Te 

(2c) 

(rr  =  Ci7K-'  +  C27-c-l  +  <pnATe 

(2d) 

where  the  radial  coordinate  was  normalized  by  the  inner  radius  ri9  and  denoted  by  r  .  The 
symbols  k  and  (p\  are  defined  in  terms  of  the  material  properties: 


K  - 


(p\  —  kS00  4*  S6r ,  (p2  —  fcS&0  S0r , 


(3a) 


<Pt\  ~  «  Q,  (pT1  —  (Srr  +  ^$r)<PT\  + 

^66  ^  rr 


(3b) 


Unknown  constants  C\  and  C 2  are  determined  from  the  boundary  conditions:  zero  radial 
stresses  at  the  inner  and  outer  surfaces  in  Eq.  (2d).  For  a  cylinder  of  multi  rings,  the  radial 
stress  and  displacement  continuity  conditions  should  be  met  at  each  ring  interface.  It  is 
noted  that  every  radial  coordinate  in  Eq.  (2)  is  expressed  in  terms  of  the  normalized 
radius  F .  We  can  thus  conclude  that  the  residual  stresses  are  the  same  for  the  cylinders 
with  the  geometric  similarity. 
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3.  EXPERIMENTS 


Three  different  winding  angles,  i.e.,  ^=90  (hoop  winding),  75  and  60  degrees  were 
considered  in  this  study  in  order  to  investigate  the  effects  of  the  winding  angles  on  the 
residual  strains  and  through-thickness  properties  of  the  cylinders.  The  cylinders  are 
referred  to  as  H90,  H75,  and  H60,  respectively.  The  ratio  of  the  outer  radius  to  the  inner 
radius  of  all  cylinders  is  given  as  3:  they  have  the  same  inner  radius  of  37.5  mm  and  outer 
radius  of  112.5  mm.  The  cylinder  specifications  were  summarized  in  Table  1. 

Manufacturing  Thick  Helical  Wound  Cylinders 

E-Glass  fibers  (RS2300-366,  Owens  Coming)  were  wound  over  the  steel  mandrel 
of  37  mm  outer  radius,  following  wet  impregnation  through  a  resin  bath.  Resin  system 
was  prepared  by  mixing  epoxy,  catalysis  and  accelerator  with  the  mass  fraction  ratio  of 
60:35:5.  The  resin  viscosity  was  maintained  the  lowest  possible  so  as  to  minimize  any 
possible  air  bubbles  that  cause  voids.  During  the  whole  filament  winding  process,  the 
fiber  tension  was  kept  as  ION,  and  the  revolution  speed  of  the  mandrel  as  50  rpm.  The 
nose  which  fed  the  fibers  was  axially  traveled  at  the  speed  regulated  to  generate  the 
desired  helical  angles.  Filament  winding  process  was  then  followed  by  curing  in  the 
autoclave  with  the  cure  cycle  as  shown  in  Figure  2.  In  order  to  ensure  the  axisymmetric 
material  properties,  the  one  end  of  the  horizontal  mandrel  shaft  was  connected  to  the 
motor  and  revolved  with  the  rotational  speed  of  15  rpm  during  the  whole  curing  process. 

Measurement  of  Through-Thickness  Mechanical  Properties 

After  cure  and  subsequent  removal  from  the  mandrel,  the  cylinder  was  sliced  and 
machined,  yielding  the  length  of  25.4  mm.  Cylinders  were  then  cut  with  a  water  cooled 
diamond  saw  into  2.5  mm  thick  specimens  for  radial  tensile  tests  as  shown  in  Figure  3. 
Two  strain  gages  were  bonded  onto  both  surfaces  and  tensile  loads  were  applied  using 
Tinius  Olsen  T1000  universal  tester  with  1  mm/min  displacement  rate. 

Specimens  for  measuring  the  CTE  in  the  radial  and  circumferential  directions  were  also 
prepared  from  the  cylinders  as  shown  in  Figure  4.  For  each  specimen,  two  high 
temperature  strain  gages  (WK-06-062AP-350,  Micro-Measurement)  were  bonded  onto 
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temperature, 


Figure  2.  Temperature  cure  cycles. 


Figure  3.  Specimen  configuration  for  tensile  testing. 


9 


(a)  radial  direction  (b)  hoop  direction 

Figure  4.  Pie  shaped  Specimen  to  measure  the  CTEs  of  Glass/Epoxy  cylinder. 


Figure  5.  Three  specimens  H90,  H75  and  H60  before  and  after  split  tests. 
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both  surfaces  of  the  specimen.  For  a  full  bridge  circuit,  the  other  two  strains  were 
bonded  to  an  ultra-low  expansion  titanium  silicate  for  temperature  compensation  of 
dummy  gages  [15].  An  adhesive  M-Bond  610  was  used  for  high  temperature  adhesive. 
The  gage  excitation  voltage  was  kept  as  2.0  V  to  avoid  self  heating  due  to  the  current. 
The  specimens  and  dummy  gages  were  heated  from  20  °C  to  80  °C .  The  averaged  slope 
of  the  measured  thermal  strains  over  the  temperature  range  was  taken  as  the  CTE. 

Measurement  of  Residual  Strains 

The  main  purpose  of  this  research  is  to  determine  the  residual  strains  in  filament 
helical  wound  cylinders  according  to  winding  angles.  Before  measuring  the  residual 
strains  by  the  destructive  split  ring  methods,  we  performed  a  non-destructive  Liquid 
Penetration  (LP)  test  to  examine  the  cracks  or  delaminations  exposed  to  the  surfaces  [19]. 
From  the  previous  experiments  [6],  the  rings  with  the  lager  radius  ratios  can  have  the 
delaminations  due  to  the  residual  strains  generated  in  the  cure  process.  This  LP  test  can 
be  easily  performed  to  visualize  any  surface  crack,  and  its  standard  procedures  are  well 
described  in  ASTM  E165  and  ISO  3879.  Once  LP  tests  are  done,  side  surfaces  were 
polished.  Two  sets  of  strain  gages,  each  for  radial  and  hoop  direction,  were  then  bonded 
on  both  surfaces  of  the  cylinder  as  shown  in  Figures  5.  Each  set  consisted  of  from  5  to  15 
stain  gages  depending  on  the  size  of  the  cylinders.  The  overall  residual  strains  were  first 
released  by  splitting  the  cylinder  in  the  radial  direction  as  shown  in  Figure  5. 

4.  RESULTS 

Young’s  Moduli,  tensile  strengths,  and  CTE  in  the  through-thickness  direction 
(radial  direction)  of  the  helically  wound  structure  for  90,  75,  and  60  degree  helical  angles 
were  measured  and  listed  in  Table  3,  where  the  average  values  of  five  specimens  with  the 
standard  deviation  were  shown.  The  through-thickness  properties  are  slightly  changed  as 
the  winding  angles  decrease:  the  radial  Young’s  Moduli  for  90,  75,  and  60  winding  angles 
were  measured  as  10.41,  10.56,  and  10.82  GPa;  the  coefficients  of  thermal  expansion  as 
28.08,  29.62,  and  30.04  X  10'6/°C;  and  the  tensile  strength  as  29.14,  25.59,  and  24.70  MPa. 
Such  decrease  along  with  the  winding  angles  are  thought  to  be  mostly  due  to  the 
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interlaminar  shear  stresses  at  the  free  edges  generated  by  the  different  contractions  along 
the  transverse  fiber  and  matrix  directions.  Another  source  may  be  from  the  matrix  rich 
zones  at  the  interface  between  layers  of  the  angle-ply  laminates.  A  further  investigation 
should  be  performed  to  understand  such  variations  using  detailed  analysis  such  as 
micromechanics. 

One  of  major  objectives  in  the  current  research  is  to  measure  and  predict  the 
residual  strains  within  helically  wound  structures.  The  residual  radial  and  hoop  strains 
measured  from  the  cylinders  of  three  different  winding  angles  90,  75,  and  60  degree  are 
shown  as  square  dots  in  Figure  6.  The  data  show  the  measured  strains  after  split  minus 
those  before  split  and  thus  represent  the  residual  strains  present  in  the  structures.  Whereas 
the  residual  strains  in  the  circumferential  directions  are  slightly  increasing  with  the 
decrease  of  winding  angles,  the  residual  strains  in  the  radial  direction  are  significantly 
decreasing,  i.e.,  1.23,  0.99,  and  0.54  X  10"3  for  90,  75,  and  60  winding  angles.  The 
maximum  radial  and  hoop  strains  versus  the  winding  angles  are  shown  in  Figure  7.  The 
decrease  in  the  radial  strains  can  be  very  important  particularly  in  a  flywheel  rotor  since 
the  tensile  radial  stresses  are  further  developed  due  to  the  centrifugal  forces.  An  effective 
temperature  change  ATe=-140°C  is  used  in  the  two-dimensional  analysis  and  three 
dimensional  FEM.  The  material  properties  of  E-Glass/Epoxy  used  in  the  analysis  are 
listed  in  Table  2.  The  results  are  in  good  agreement  with  the  experimental  data.  The 
discrepancies  between  the  two-dimensional  analysis  and  the  FEM  increase  as  the  winding 
angle  decrease.  It  is  probably  because  the  z-directional  effects  neglected  in  the  two- 
dimensional  analysis  emerge  in  the  three-dimensional  FEM.  The  variation  of  the 
maximum  radial  strains  along  with  the  outer  radius  of  the  cylinders  is  obtained  for  three 
winding  angles  and  shown  in  Figure  8.  It  is  observed  that  the  maximum  radial  strains  are 
significantly  affected  by  the  outer  radius.  As  noticed  in  the  analysis,  the  strains  are 
dependent  upon  only  the  radius  ratio,  which  means  that  the  residual  strains  are  the  same 
for  the  cylinders  with  the  geometric  similarity. 
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(a)  radial  and  hoop  strains  of  specimen  H90 


(b)  radial  and  hoop  strains  of  specimen  H75 


(c)  radial  and  hoop  strains  of  specimen  H60 


Figure  6.  Distribution  of  radial  and  hoop  strains  of  helically  wound  cylinders  H90,  H75, 
and  H60  along  the  radius:  the  experimental  data  were  compared  with  the  two- 
dimensional  analysis  and  three-dimensional  FEM. 
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winding  angle,  degree 
(a)  radial  strains 

Figure  7 .  Variation  of  the  maximum  radial 
angles. 


winding  angle,  degree 

(b)  hoop  strains 

and  hoop  residual  strains  with  the  winding 


outer  radius,  rQ  outer  radius,  r0 

(a)  max.  radial  strains  (b)  max.  hoop  strains 


Figure  8.  Variation  of  maximum  radial  and  hoop  residual  strains  along  the  outer  radius 
ratios  for  three  different  winding  angles  90,  75,  and  60  degrees. 
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radial  stress,  MPa  radial  stress,  MPa  radial  stress,  MPa 


(a)  radial  and  hoop  stresses  of  specimen  H90 


(b)  radial  and  hoop  stresses  of  specimen  H75 


radial  coordinate,  r  radial  coordinate,  y 


(c)  radial  and  hoop  stresses  of  specimen  H60 

Figure  9.  Distribution  of  radial  and  hoop  stresses  of  the  helically  wound  rings  along  the 
radius  calculated  by  the  present  analysis  and  FEM. 
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The  distribution  of  the  radial  and  circumferential  stresses  were  also  obtained  and 
shown  in  Figure  9.  The  maximum  circumferential  stresses  always  occur  at  the  inner 
surfaces  whereas  the  maximum  radial  stresses  occur  inside  the  cylinder.  The  maximum 
radial  stress  of  the  H90  cylinder  reaches  as  high  as  40%  of  the  radial  strength.  The 
circumferential  stresses,  tensile  at  the  inner  surface  and  compressive  at  the  outer  surface, 
generate  the  bending  moment  inside  the  radial  surface,  which  tends  to  rotate  the  radial 
surface  once  split.  Such  bending  moment  is  also  related  to  the  tensile  radial  stresses 
noticed  in  this  study. 


5.  CONCLUSION 

We  have  measured  the  through-thickness  material  properties  and  the  residual  strains 
of  helically  thick  wound  cylinders.  Whereas  the  radial  stiffness  remains  almost  the  same 
from  90  to  60  degree  winding  angles,  about  7  %  increase  of  CTE  and  about  17  % 
decrease  of  the  radial  strengths  was  observed.  Such  variations  of  the  material  properties 
with  the  winding  angles  are  believed  to  be  due  to  interlaminar  shear  stresses  at  the  free 
edges  generated  by  the  different  contractions  along  the  transverse  fiber  and  matrix 
directions  in  the  test  specimens.  It  is  believed  that  such  variations  can  be  better  explained 
using  a  detailed  analysis  such  as  micromechanics  approach.  The  maximum  residual 
strains  in  the  radial  direction  were  observed  to  decrease  by  about  54%  from  H90  to  H60 
cylinders  whereas  the  maximum  residual  strains  in  the  hoop  increase  by  about  25%.  Such 
decrease  in  the  radial  direction  can  be  effectively  utilized  in  the  application  for  flywheel 
rotors  in  which  the  radial  strengths  are  crucial  in  supporting  the  centrifugal  forces.  Design 
of  such  rotors  often  requires  an  optimization  procedure  together  with  the  structural 
analysis  that  should  consider  the  residual  strains. 
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Table  1.  Specification  of  the  helically  wound  rings  for  measurement  of  residual  strains. 


Name 

Specimen  H90 

Specimen  H75 

Specimen  H60 

Ratio  (ro/n) 

3.0 

3.0 

3.0 

r\  (mm) 

37.5 

37.5 

37.5 

r0  (mm) 

110.5 

110.5 

110.5 

Material 

E-Glass/Epoxy 

E-Glass/Epoxy 

E-Glass/Epoxy 

Length  (mm) 

25.4 

25.4 

25.4 

Winding  angle 

90°  (hoop) 

75° 

60° 

Table  2.  Material  properties  of  E-Glass/Epoxy  used  in  the  analysis 


Material  Properties 

Symbols 

Unit 

E-Glass/Epoxy 

Longitudinal  Young’s  Modulus 

Ex 

GPa 

38.6 

Transverse  Young’s  Modulus 

Ey=  Er 

GPa 

See  Table  3. 

Stiffness  Ratio 

K 

2.16 

Longitudinal  Poisson’s  Ratio 

Vx 

0.26 

Longitudinal  Shear  Modulus 

Gx 

GPa 

4.14 

Longitudinal  CTE 

ax 

10'6/°C 

6.35 

Transverse  CTE 

ay=ar 

10'6/°C 

See  Table  3. 

Longitudinal  Tensile  Strength 

X 

MPa 

1062 

Longitudinal  Comp.  Strength 

X' 

MPa 

610 

Transverse  Tensile  Strength 

Y 

MPa 

See  Table  3. 

Transverse  Comp.  Strength 

Y’ 

MPa 

118 

Longitudinal  Shear  Strength 

S 

MPa 

72 

Table  3.  Through-thickness  material  properties  measured  from  helically  wound  E- 

Glass/Epoxy  rings. 


Winding 

E-Glass/Epoxy 

Angle 

Er,  GPa 

Y,  MPa 

ccr,  10'6/°C 

90°(hoop) 

10.4110.67 

29.1411.87 

28.0810.31 

75° 

10.56+1.13 

25.5913.44 

29.6210.63 

60° 

10.8212.82 

24.7014.32 

30.0410.22 

Average 

10.59 

26.47 

29.24 
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PART  II.  Effects  of  Rotor  Sizes  and  Epoxy  System  on  the  Process- 
Induced  Residual  Strains  within  Multi-ring  Composite  Rotors 

1.  INTRODUCTION 

The  fiber  reinforced  composite  materials  have  been  successfully  used  for  the 
flywheel  rotors.  It  is  partly  due  to  the  anisotropic  properties  which  results  in  the  reduction 
of  the  stresses  generated  by  the  rotational  centrifugal  forces  [1-4].  Different  fibers  such  as 
Glass  and  Graphite  fibers  can  be  further  combined  into  multi-ring  rotors  to  reduce  the 
stresses  and  increase  even  more  performance  of  the  flywheel  rotors  [5~7].  On  the  other 
hand,  however,  the  anisotropic  stiffness  and  coefficients  of  thermal  expansion  differing  in 
the  radial  and  circumferential  directions  inevitably  generate  the  residual  tensile  radial 
strains,  which  significantly  deteriorate  the  performance  of  the  rotor  since  the  radial 
tensile  stresses  are  mostly  sustained  by  the  relatively  weak  matrix.  The  process-induced 
residual  strains  can  be  so  high  that  it  leads  to  delamination  within  thick  filament  wound 
rotors  [8].  Those  residual  strains  are  generally  proportional  to  the  temperature  drop  and 
chemical  shrinkage  experienced  during  the  curing  process  [9~14],  Thus  the  temperature 
profile  in  a  curing  cycle  is  a  significant  process  parameter  affecting  the  strength  and 
integrity  of  the  flywheel  rotor.  Understanding  the  effects  of  the  resin  systems  on  the 
residual  strains  within  thick  multi-ring  rotors  is  critical  in  design  of  the  composite 
flywheel  rotors  for  larger  energy  storage  capacities. 

Prediction  of  the  residual  strains  generated  during  the  curing  process  in  thick 
composite  rotors  generally  requires  an  understanding  of  very  complicated  multi-physics 
such  as  a  detailed  universal  rheological  model  including  cure  shrinkage  as  well  as 
thermal  strains,  the  process  of  gelation,  cure  kinetics,  the  exotherm  and  interaction 
between  the  rotor  and  the  mandrel.  Several  investigators  have  developed  theoretical 
models  of  the  filament  winding.  Among  several  attempts,  Springer  et  al.  modeled  the 
processes  applicable  to  the  moderately  thick  rings  by  establishing  submodels,  each  for 
thermochemical  behavior,  fiber  motion,  void,  and  stress,  and  developed  a  cure  simulation 
computer  program  called  WINDTHICK  [11-13].  A  simplified  approach  using  an 


20 


effective  temperature  drop  which  was  verified  using  the  WINTHICK  was  successfully 
used  to  predict  the  redial  strains  within  rotors  of  Graphite/Epoxy  [8].  However,  the 
effects  of  cure  cycles  and  Epoxy  system  on  the  residual  strains  have  never  been 
pertinently  addressed. 

Multi-ring  type  material  hybridization  in  flywheel  rotors  has  been  very  effectively 
used  to  reduce  the  tensile  stresses  in  the  radial  direction  generated  by  centrifugal  forces 
[5—7] .  Generally,  arrangement  of  softer  materials  at  inner  sides  and  harder  materials  at 
outer  sides  results  in  the  reduction  of  radial  tensile  stresses.  Several  rotors  were  optimally 
designed  which  maximize  the  overall  performance  of  the  flywheel  energy  storage  system. 
A  structural  analysis  in  the  design  procedure  considered  both  centrifugal  forces  and  the 
temperature  drop  that  occurred  during  the  curing  process  [6,  7].  However,  the  residual 
stresses  within  thick  multi-ring  rotors  have  never  been  experimentally  measured  and 
quantified  so  far. 

In  this  study,  the  effects  of  the  different  epoxy  systems  cured  at  different 
temperature  profiles  and  the  dimensional  size  on  the  residual  strains  within  the  hybrid 
multi-ring  rotors  will  be  discussed.  Rotors  of  various  sizes  will  be  manufactured  using 
Glass  and  Graphite  fibers  with  different  epoxy  systems.  The  radial  stiffness  and  strength 
and  the  coefficients  of  thermal  expansion  will  be  measured  to  evaluate  the  characteristics 
of  the  epoxy  systems.  Anon-destructive  liquid  penetration  tests  [15]  will  be  performed  to 
visually  examine  the  radial  delamination,  followed  by  a  destructive  split-ring  test  to 
measure  the  residual  strains  [8,16,17].  An  effective  temperature  will  be  then  determined 
and  the  residual  stresses  and  strains  will  be  estimated  using  a  stress  analysis.  The  main 
purpose  of  this  research  is  to  find  an  epoxy  system  which  will  lower  residual  strains 
within  the  hybrid  rotors  while  maintaining  the  radial  strengths  as  high  as  the  conventional 
high  temperature  cured  rotors. 

2.  STRESS  AND  STRENGTH  ANALYSIS 

The  radial  delamination  can  occur  within  a  ring  due  to  residual  stresses.  A 
maximum  stress  criterion  is  used  for  the  creation  of  delamination  [19]: 
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(1) 


Delamination  occurs  when  the  delamination  index  Da  reaches  the  value  of  one.  Once  a 
delamination  occurs,  the  ring  is  assumed  be  split  into  two  rings.  The  radial  location 
ra  max ,  where  the  maximum  radial  stress  occurs,  can  be  determined  from  the  stress 

analysis  using  Eqs  (2~3)  in  Part  I  and  the  boundary  conditions  for  hybrid  rotors.  The 
maximum  temperature  drops  ATD  which  cause  an  initial  delamination  can  be  also 
obtained  using  Eqs.  (1—4)  in  Part  I.  The  hybrid  multi -ring  rotors  of  the  geometric 
similarity  experience  a  delamination  at  the  same  effective  temperature  drop. 

3.  EXPERIMENTS 

Measurement  of  Mechanical  Properties 

Three  different  matrix  systems  were  considered  in  this  study  in  order  to  investigate 
the  effects  of  the  matrix  systems  on  the  residual  strains  within  the  rotors.  Their  chemical 
compound  names  and  mixing  ratios  were  listed  in  Table  1,  each  denoted  by  high,  mid, 
and  low  temperature  epoxy  systems  (HTES,  MTES,  LTES,  respectively),  and  the 
corresponding  cure  cycles  recommended  by  the  resin  suppliers  were  shown  in  Figure  1. 
The  LTES  can  be  cured  alternatively  at  room  temperature  for  about  48  hours. 

Specimens  were  prepared  for  measurement  of  the  mechanical  properties  of  the  pure 
resin  systems.  First,  steel  mold  was  machined  and  specimen  with  each  epoxy  system  was 
cured  according  to  the  given  cure  cycles.  Strain  gages  (AE-XX-S30N-350,  CAS)  were 
then  attached  to  the  specimens,  and  Young’s  Moduli  and  tensile  strengths  of  the  pure 
matrix  systems  were  measured  according  to  the  JIS  K  6850  specification. 

To  evaluate  the  characteristics  of  the  epoxy  systems,  the  radial  stiffness  and 
strength  and  the  coefficients  of  thermal  expansion  were  measured  from  the  Glass  and 
Graphite  rotors.  The  steel  mandrel  of  37.5  mm  outer  radius  was  first  machined  with  two 
side  circular  plates  of  the  outer  radius  of  150  mm  and  the  thickness  of  10  mm  attached  for 
guiding  the  winding  fiber  to  produce  the  rotor  of  length  of  200  mm.  Graphite  (T300-12K, 
TORAY)  fiber  and  E-Glass  fiber  (RS2300-366,  Owens  Coming)  were  used,  and  the  fiber 
tension  of  8-10  N  was  applied  during  the  whole  filament  winding  process.  The  mandrel 


22 


was  rotating  with  the  revolution  speed  of  50  rpm,  and  the  nose  which  fed  the  fibers 
axially  traveled  at  the  speed  of  1  m/min.  The  wet  filament  winding  process  was  then 
followed  by  a  curing  process  in  the  autoclave  with  the  temperature  cycles  as  shown  in 
Figure  1.  In  order  to  ensure  the  axisymmetric  material  properties,  the  one  end  of  the 
horizontal  mandrel  shaft  was  connected  to  the  motor  and  revolved  with  the  rotational 
speed  of  15  rpm  during  the  whole  curing  process.  After  cure  and  subsequent  removal 
from  the  mandrel,  the  rotors  were  sliced  and  machined  using  a  water  cooled  diamond  saw 
into  15  mm  length  rotor.  Rotors  were  then  cut  with  a  water  cooled  diamond  saw  into  2.5 
mm  thick  specimens  for  radial  tensile  tests  as  shown  in  Figure  2.  The  two  strain  gages 
were  bonded  onto  both  surfaces  to  measure  the  radial  strains  in  the  tension  tests. 

Specimens  for  measuring  the  CTE  in  the  radial  and  circumferential  directions  were 
also  prepared  from  the  rotors  as  shown  in  Figure  3.  For  each  specimen,  two  high 
temperature  strain  gages  (WK-06-062AP-350,  Micro-Measurement)  were  bonded  onto 
both  surfaces  of  the  specimen.  For  a  full  bridge  circuit,  the  other  two  strains  were  bonded 
to  an  ultra-low  expansion  titanium  silicate  for  temperature  compensation  of  dummy 
gages  [21].  An  adhesive  M-Bond  610  was  used  for  high  temperature  adhesive.  The  gage 
excitation  voltage  was  kept  as  2.0  V  to  avoid  self  heating  due  to  the  current.  The 
specimens  and  dummy  gages  were  heated  from  20  °C  to  80  °C .  The  averaged  slope  of  the 
measured  thermal  strains  over  the  temperature  range  was  taken  as  the  CTE. 

Measurement  of  Residual  Strains 

The  main  purpose  of  this  research  is  to  determine  the  residual  strains  in  filament 
wound  hybrid  rotors  according  to  the  matrix  systems  and  the  rotor  sizes.  Three  two-ring 
rotors  of  three  different  sizes  (outer  radius  ratios  70  =1.5,  2.0  and  2.5)  for  each  epoxy 
system  (HTES,  MTES,  LTES)  were  prepared  using  the  wet  filament  winding  methods. 
The  same  mandrel,  side  plates,  and  winding  methods  were  used  as  for  the  rotors  for 
measurement  of  the  material  properties:  all  rotors  had  the  inner  radii  of  37.5  mm.  All 
rotors  were  hoop  wound  (90  degree  winding  angle)  using  E-Glass  fibers  for  the  inner  ring 
and  T300-12K  graphite  fiber  for  the  outer  ring  with  the  same  thickness.  The  rotors  of 
different  sizes  were  denoted  respectively  by  R1.5,  R2.0,  and  R2.5,  and  the  specifications 
were  summarized  in  Table  5.  Each  series  of  rotors  were  then  cured  in  an  autoclave 
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according  to  the  epoxy  systems,  and  Rx.x-H,  M  or  L  are  used  to  denote  a  rotor  of  radius 
ratio  x.x  with  high,  middle  or  low  temperature  epoxy  systems,  respectively.  The  multi¬ 
ring  rotors  were  wet-wound  and  cured  in  the  same  way  previously  described.  The  rotors 
were  then  removed  from  the  mandrel,  and  machined  into  25.4  mm  length  rotor. 

Before  measuring  the  residual  strains,  we  performed  a  non-destructive  Liquid 
Penetration  (LP)  test  to  examine  the  delaminations  exposed  to  the  surfaces  [15].  From  the 
previous  experiments  [8],  it  was  observed  that  the  rotors  with  the  large  radius  ratios  can 
have  the  delaminations  due  to  the  residual  strains  generated  in  the  cure  process.  This  LP 
test  can  be  easily  performed  to  visualize  any  surface  crack,  and  its  standard  procedures 
are  well  described  in  ASTM  El 65  and  ISO  3879.  To  be  brief,  the  surfaces  are  first 
cleaned;  a  colored  liquid  is  applied  to  the  surface  and  it  may  penetrate  into  the  cracks  by 
capillary  action;  after  the  removal  of  any  excess,  an  absorption  white  layer  is  applied, 
drawing  the  colored  liquid  out  of  the  cracks  and  so  making  them  visible.  The  penetrant, 
removal  and  developer  were  supplied  from  NAWOO  Tech  Research  Ltd.. 

Once  LP  tests  are  done,  the  surfaces  were  cleaned  using  sand  papers,  and  two  sets 
of  strain  gages,  each  for  radial  and  hoop  direction,  were  bonded  on  both  surfaces  of  the 
rotor  as  shown  in  Figures  4.  Each  set  consisted  of  from  5  to  15  stain  gages  depending  on 
the  size  of  the  rotors.  The  overall  residual  strains  within  the  rotors  were  released  by 
splitting  the  ring  in  the  radial  direction  as  shown  in  Figure  4.  A  radial  slice  of  about  10 
degree  containing  each  row  of  strain  gages  was  then  cut  to  completely  release  the  residual 
strains  for  hybrid  rotors.  The  second  cuttings  were  necessary  since  the  first  cut  may  not 
completely  release  the  residual  strains  in  the  hybrid  multi-ring  rotors.  The  readouts  before 
the  first  cutting  were  set  to  zero,  and  then  recorded  after  first  and  second  cutting. 
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Figure  1.  Three  different  cure  cycles  respectively  for  HTES,  MTES, 


and  I 


Figure  2.  Specimen  configurations  for  tensile  tests. 
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(a)  radial  direction 


(b)  hoop  direction 


Figure  3.  Pie  shaped  specimen  with  the  stain  gages  for  measurement  of  the  CTEs. 


(a)  specimen  before  split 


(b)  specimen  after  1st  split  (c)  specimen  after  2nd  split 


Figure  4.  The  R2.5-H  rotor  before  and  after  split  in  the  split-ring  test. 


4.  RESULTS 


Material  Properties  of  Three  Resin  Systems 

Young’s  Moduli  and  tensile  strengths  of  the  specimens  of  pure  matrix  systems  were 
measured  and  their  average  values  with  the  standard  deviation  are  listed  in  Table  3.  It  is 
noted  that  the  pure  MTES  has  the  highest  tensile  strength  (60  MPa)  while  the  pure  HTES 
has  the  highest  Young’s  Modulus  (4.88  GPa).  Exact  reasons  are  unknown  at  this  point, 
but  this  MTES  can  be  an  excellent  candidate  for  the  rotor  since  lower  cure  temperature 
can  generally  lower  the  residual  strains  within  the  rotors. 

The  resin  system  can  directly  affect  the  radial  Young’s  modulus  and  tensile  strength 
of  the  rotors,  and  their  properties  were  measured  from  specimen  cut  from  the  rotors  with 
three  different  epoxy  systems.  Those  measured  values  are  listed  in  Table  3.  Among  the 
epoxy  systems,  the  HTES  yields  the  highest  stiffness  in  both  rotors  of  Glass  and  Graphite 
fibers:  about  5%  higher  than  the  rotors  with  MTES.  However,  the  MTES  rotors  shows 
about  2~7%  higher  strength  than  the  HTES  rotors.  The  LTES  in  both  rotors  of  Glass  and 
Graphite  exhibits  about  35%  lower  strength  than  other  epoxy  systems.  However,  the 
Young’s  Moduli  shows  little  variation  over  the  considered  epoxy  systems. 

Besides  to  the  stiffness  anisotropy,  the  difference  between  the  CTE  in  the  radial  and 
circumferential  directions  is  causing  the  residual  strains  as  noticed  in  Eqs.  (2)  and  (3)  in 
Part  I.  The  measured  CTE  of  rotors  of  Glass  and  Graphite  fibers  were  in  Table  4.  As 
expected,  the  near  zero  thermal  expansion  in  the  circumferential  direction  was  observed 
in  the  Graphite  rotors.  In  the  case  of  rotors  of  Glass  fiber,  however,  the  CTE  in  the 
circumferential  direction  is  measured  as  6.35  x  1  (T6  /  °C ,  which  is  about  20%  of  the  CTE 
in  the  radial  direction.  These  CTE  were  measured  from  the  rotors  with  the  MTES,  and 
assumed  to  be  the  same  for  other  matrix  systems. 

Effects  of  Rotor  Sizes  on  the  Residual  Strains 

One  of  the  main  purposes  of  the  current  research  is  to  investigate  the  effects  of 
rotor  sizes  of  the  multi-ring  rotors  on  the  residual  strains.  The  residual  radial  and  hoop 
strains  measured  from  the  rotors  of  three  different  sizes,  i.e.,  R1.5-H,  R2.0-H,  and  R2.5-H 
(all  with  HTES)  are  shown  in  Figure  5.  The  data  show  the  strains  measured  before  split 
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minus  those  before  split  and  thus  represent  the  residual  strains  present  in  the  cured  rotors. 
Changes  in  the  radial  and  circumferential  strain  distributions  are  due  to  the  material 
change  from  E-Glass  to  Graphite  fibers.  Differently  from  the  rotor  R1 .5-H,  delaminations 
are  noticed  in  rotors  R2.0-H  and  R2.5-H  from  the  LP  tests  as  shown  in  Figure  6.  On 
splitting  the  specimen,  the  delaminated  surfaces  were  thus  contracted  in  the  radial 
direction,  and  discontinuously  deformed  in  the  hoop  direction,  as  shown  in  Figures  5(a) 
and  (b).  The  effective  temperature  A Te  =-140°C  was  first  determined  from  the  rotor 

R1.5-H,  and  applied  to  the  other  rotors  R2.0-H  and  R2.5-H.  The  calculated  residual 
strains  are  shown  in  Figures  5(a)  and  (b).  The  material  properties  of  E-Glass/Epoxy  and 
T300/Epoxy  used  in  the  analysis  are  listed  in  Table  1  and  2  [19].  The  measured  properties 
were  used  for  the  epoxy  dominant  material  properties.  The  calculated  strains  are  in  good 
agreement  with  the  measured  residual  strains. 

The  residual  stresses  of  the  HTES  rotors  with  ATe  =  -140°C  are  calculated  and 
shown  in  Figure  7.  The  effective  temperature  drops  causing  the  delamination  in  each 
specimen  are  calculated  using  Eq.  (1):  ATD  =  -263,  -138  and  -105  °C  respectively  for 
specimens  R1.5-H,  R2.0-H  and  R2.5-H.  Notice  that  |A7^|  =  140  >  |ATfl|  for  specimens 
R2.0-H  and  R2.5-H,  indicating  occurrence  of  delaminations  (  Da  >1).  Actually,  the 
effective  temperature  drop,  A Te  =-140 °C  yields  the  damage  indices  Da  of  0.53,  1.01 
and  1 .32  respectively  for  R1 .5-H,  R2.0-H  and  R2.5-H  rotors. 

Effects  of  Resin  Systems  on  the  Residual  Strains 

In  order  to  reduce  the  residual  strains,  we  studied  two  other  epoxy  systems  which 
can  be  cured  at  lower  temperature  than  the  HTES.  Rotors  were  manufactured  accordingly, 
and  the  residual  strains  were  measured  as  described  in  the  previous  section  and  shown  in 
Figures  8  for  three  different  rotor  sizes.  Circles  represent  the  data  for  the  HTES,  squares 
for  the  MTES,  and  triangles  for  the  LTES;  the  data  points  are  average  values  of  least  four 
experimental  data. 

It  is  noted  in  the  R1.5  rotors  that  the  residual  strains  with  the  MTES  were  reduced 
by  about  65  %  of  those  with  HTES  as  shown  in  Figure  8(a).  Those  reductions  in  the 
residual  strains  were  obviously  resulted  from  the  lower  temperature  profile  of  MTES  cure 
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cycle  than  HTES.  The  residual  strains  of  rotors  with  LTES  are  as  low  as  those  with 
MTES,  but  the  radial  strengths  attained  only  60%  of  those  with  HTES.  Effective 
temperature  changes  were  determined  as  A 7’  =  — 50°C  for  rotor  R1.5-M,  and 

ATe  =  -35  °C  for  rotor  R1.5-L. 

The  measured  residual  strains  of  rotors  R2.0  and  R2.5  were  shown  respectively  in 
Figure  8(b)  and  (c).  Whereas  we  have  observed  the  delaminations  in  rotors  of  R2.0-H  and 
R2.5-H,  no  delaminations  were  detected  in  rotors  R2.0-M,  L  and  R2.5-M,  L.  The 
effective  temperatures  determined  from  rotors  of  R  1.5  were  used  to  calculate  the  residual 
strains  for  rotors  of  bigger  sizes  and  shown  in  Figures  8(b)  and  (c).  The  calculated  values 
are  in  good  agreement  with  the  measured  data.  It  can  be  concluded  that  an  effective 
temperature  can  be  used  to  predict  the  residual  strains  of  the  rotors  for  different  sizes  as 
long  as  they  are  made  of  the  same  material  systems.  This  observation  agreed  with  the 
derived  results  in  Section  2,  where  the  residual  strains  and  stresses  were  expressed  in 
terms  of  dimensionless  radial  coordinate. 

The  residual  stresses  corresponding  to  the  residual  strains  were  calculated  and 
shown  in  Figure  9.  For  rotors  of  R2.0-H  and  R2.5-H  which  experienced  delaminations, 
the  stresses  without  consideration  of  delaminations  were  shown  to  compare  the  maximum 
possible  stresses.  As  in  the  residual  strains,  significant  stress  reductions  were  noticed  in 
rotors  with  MTES  and  LTES  for  all  sizes  of  rotors.  The  discontinuities  in  the 
circumferential  stress  distribution  are  obviously  due  to  the  material  change  from  Glass  to 
Graphite  fibers.  It  is  noted  that  the  maximum  radial  stresses  always  occur  at  the  interface 
between  the  two  rings.  This  observation  can  be  predicted  using  the  analysis  presented  in 
Section  2. 

The  locations  of  cr.  max  versus  the  radius  ratio  Fa  are  calculated  and  shown  as  solid 

lines  in  Figure  10(a).  Since  the  maximum  radial  stresses  occurs  at  the  ring  interface  and 
each  ring  has  the  same  thickness,  the  radial  locations  coincide  with  the  centerline  of  the 
overall  rotor  thickness.  The  effective  temperature  changes  which  cause  the  delaminations, 
denoted  by  ATD  (allowable  effective  temperature  drops  before  delamination),  are 
calculated  in  terms  of  the  radius  ratios  rB  for  each  epoxy  system  and  shown  in  Figure 
10(b);  ATd  =  -263,  -138  and  -105  °C  respectively  for  rotors  R1.5-H,  R2.0-H  and  R2.5-H, 
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ATD-  -271,  -143  and  -110°C  for  rotors  R1.5-M,  R2.0-M  and  R2.5-M,  and  ATD  =  -181,  - 
96  and  -73  °C  for  rotors  R1.5-L,  R2.0-L  and  R2.5-L.  The  A TD  for  LTES  is  quite  lower 
than  others  since  their  lower  radial  strengths.  Remember  the  ATe ’s  for  HTES,  MTES  and 

LETS  were  determined  respectively  as  -150,  -50  and  -35  °C.  Similarly,  the  variation  of 
the  delaminations  indices  D ^  for  each  epoxy  system  over  the  outer  radii  r0  are 

calculated  and  shown  in  Figure  10(c).  It  is  once  more  observed  that  the  MTES  has  the 
lowest  delamination  indices  due  to  the  facts  that  the  MTES  rotors  have  lower  temperature 
drop  than  the  HTES  rotors  and  have  the  higher  radial  strength  than  the  LTES  rotors.  In 
other  words,  the  higher  radial  strength  together  with  the  lower  residual  stresses  can  be 
achieved  in  the  MTES  rotors. 


5.  CONCLUSION 

We  have  measured  and  predicted  the  process  induced  residual  strains  within  thick 
hybrid  rotors  for  different  matrix  systems.  As  noticed  in  the  results,  delaminations  can  be 
easily  developed  within  the  rotor  with  the  HTES  due  to  the  high  temperature  cure  cycles. 
We  have  observed  delaminations  in  rotors  of  the  radius  ratios  greater  than  2.0  when  using 
the  HTES.  The  delaminations  were  clearly  observed  using  both  liquid  penetration  and 
split  ring  tests.  The  residual  strains  can  be  reduced  significantly  using  the  LTES,  but  the 
radial  strengths  were  decreased  at  the  same  time.  Among  the  considered  epoxy  systems, 
the  MTES  which  needs  lower  cure  temperature  than  the  HTES  was  found  to  be  the  most 
effective  in  reducing  the  residual  strains  while  maintaining  as  high  radial  strength  as  the 
HTES.  The  effective  temperature  drops  in  calculation  of  the  residual  strains  was 
determined  as  -50  °C  for  MTES  rotors,  and  no  delamination  were  detected  in  rotors  of 
radius  ratios  up  to  2.5.  The  lower  residual  strains  can  increase  the  performance  of  the 
rotors:  the  damage  or  delamination  index  0.34  in  the  R2.0-M  means  that  the  residual 
strengths,  i.e.,  0.66  can  be  used  to  support  the  centrifugal  forces  when  used,  for  instance, 
as  the  flywheel  rotors. 
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radial  coordinate,  y 


radial  coordinate,  y 


radial  coordinate  is  normalized  by  the  inner  radius). 


(a)  R1.5-H  (b)  R2.0-H  (c)  R2.5-H 


Figure  6.  Liquid  penetrant  test  results  of  R1.5-H,  R2.0-H  and  R2.5-H  rotors.  The  red 
circles  indicate  the  delamination  caused  by  the  residual  stresses. 


(a)  radial  stresses  of  specimen 


(b)  hoop  stresses  of  specimen 


Figure  7.  Calculated  radial  and  circumferential  stresses  of  R1.5-H,  R2.0-H  and  R2.5-H  rotors 
without  consideration  of  delamination.  (zire=-140°C) 
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(a)  radial  and  hoop  strains  of  specimen  R1 .5 


radial  coordinate,  r 


radial  coordinate,  r 


(b)  radial  and  hoop  strains  of  specimen  R2.0 


b 


radial  coordinate,  r 


radial  coordinate,  r 


(c)  radial  and  hoop  strains  of  specimen  R2.5 


Figure  8.  Comparison  of  the  residual  strains  of  the  R1.5,  R2.0,  and  R2.5  rotors  for  three 
different  epoxy  systems  LTES,  MTES  and  HTES. 
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radial  stress,  MPa  radial  stress,  MPa  radial  stress,  MPa 


radial  coordinate,  y  radial  coordinate,  r 


(a)  radial  and  hoop  stresses  of  specimen  R1 .5 


(b)  radial  and  hoop  stresses  of  specimen  R2.0 


(c)  radial  and  hoop  stresses  of  specimen  R2.5 


Figure  9.  Calculated  residual  stresses  of  the  R1.5,  R2.0,  and  R2.5  rotors  for  three  different 
epoxy  systems  LTES,  MTES  and  HTES. 
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7.0 


outer  radius, r0 


(a)  location  of  maximum  radial  stresses 


(b)  allowable  temperature  drops  before  delamination 


outer  radius,  r0 

(c)  delamination  indices 


Figure  10.  Radial  location  of  maximum  radial  stresses,  allowable  effective  temperature  drop 
before  delamination,  and  delamination  indices  versus  the  radius  ratios  for  three  epoxy  systems 
LTES,  MTES,  and  HTES 
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Table  1.  Material  Properties  of  three  epoxy  system  considered  in  this  study. 


Epoxy  system 

Resin  Type 
(Commercial  Name) 
Mixture  Weight  Ratio 

Hardener 

(Commercial  Name) 
Mixture  Weight  Ratio 

Catalyzer 

(Commercial  Name) 
Mixture  Weight  Ratio 

Bisphenol-F 

Modified  Amine 

LTES 

(SS80-WEH) 

(SW1500TKD) 

N/A 

80% 

20% 

Bisphenol-F 

Modified  Amine 

MTES 

(SS80-WEH) 

(SS-2033H) 

N/A 

75% 

25% 

Bisphenol-F 

Acid  Anhydride 

Amine 

HTES 

(JF7006) 

(JF7006H) 

(JA001) 

60% 

35% 

5% 

Table  2.  Material  properties  of  Graphite/Epoxy  and  E-Glass/Epoxy 


Material  Properties 

Symbols 

Unit 

Graphite/Epoxy 

E-Glass/Epoxy 

Longitudinal  Young’s  Modulus 

Ee 

181 

38.6 

Transverse  Young’s  Modulus 

ii 

See  Table  3. 

Stiffness  Ratio 

K 

4.19 

2.16 

Longitudinal  Poisson’s  Ratio 

ve 

0.28 

Longitudinal  Shear  Modulus 

Gq 

GPa 

7.17 

4.14 

Longitudinal  CTE 

ae 

10'6/°C 

See  Table  4. 

Transverse  CTE 

ar 

io'Vc 

See  Table  4. 

Longitudinal  Tensile  Strength 

X 

MPa 

1500 

1062 

Longitudinal  Comp.  Strength 

X’ 

MPa 

1500 

610 

Transverse  Tensile  Strength 

Y 

MPa 

See  Table  3. 

Transverse  Comp.  Strength 

T 

MPa 

246 

118 

Longitudinal  Shear  Strength 

S 

MPa 

68.0 

72 
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Table  3.  Measured  Young’s  Moduli  and  the  tensile  strengths  of  E-Glass/Epoxy, 
Graphite/Epoxy  and  pure  epoxy  resin  specimens. 


Epoxy 

E-Glass/Epoxy 

Graphite/Epoxy 

Pure  Epoxy  Resin 

System 

Er,  GPa 

Y,  MPa 

Er,  GPa 

Y,  MPa 

Er,  GPa 

Y,  MPa 

LTES 

8.61±0.47 

19.1211.76 

12.6311.69 

4.1510.23 

49.6012.61 

MTES 

8.7610.60 

31.2713.93 

8.4410.61 

18.8813.54 

4.0110.27 

60.0314.19 

HTES 

8.8911.18 

29.1412.74 

9.0210.92 

18.4111.90 

4.8810.27 

56.7814.52 

Table  4.  Measured  CTEs  of  E-Glass/Epoxy  and  Graphite/Epoxy  specimens. 


E-Glass/Epoxy  CTE 

Graphite/Epoxy  CTE 

cte 

ar 

ae 

ar 

6.35x1 0'6/°C 

28.08xl0‘6/°C 

0.005x10‘6/°C 

25.46x10‘6/°C 

Table  5.  Specifications  of  the  hybrid  multi-ring  composite  rotors  for  measurement 
of  residual  strains. 


Name 

R1.5 

R2.0 

R2.5 

Ratio  (ro/rO 

1.5 

2.0 

2.5 

Y\  (mm) 

37.5 

37.5 

37.5 

r0(mm) 

56.25 

75 

93.75 

Material 

E-Glass/Epoxy 

+ 

E-Glass/Epoxy 

+ 

E-Glass/Epoxy 

+ 

Graphite/Epoxy 

Graphite/Epoxy 

Graphite/Epoxy 

Length  (mm) 

25.4 

25.4 

25.4 

Winding  Angle 

90°(Hoop) 

90°(Hoop) 

90°(Hoop) 

(Thickness  ratio  of  E-Glass/Epoxy  to  Graphite/Epoxy  rings  =1:1) 
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PART  III.  OPTIMIZATION  PROGRAM  FOR  FLYWHEEL 
COMPOSITE  ROTORS  -  FLYOPT 


1.  2-DIMENSIONAL  STRESS  &  STRENGTH  ANALYSIS 

For  a  numerical  calculation,  the  composite  flywheel  rotor  of  varying  material 
properties  is  divided  into  many  rings  and  each  ring  is  assumed  to  have  a  constant  material 
property  of  an  angle-ply  laminate  [(±<f>)m]  as  shown  in  Figure  1.  Based  on  the 

assumption  of  a  generalized  plane  strain  state,  all  the  rings  are  free  to  expand  or  contract 
in  the  axial  direction  but  have  the  same  axial  strain.  The  flywheel  rotor  is  also  assumed 
cylindrically  ortho  tropic  and  subject  to  axisymmetric  centrifugal  forces  due  to  rotation. 
With  the  axisymmetric  conditions,  the  circumferential  displacements  are  thus  neglected. 
All  the  displacements  and  stresses  are  thus  independent  of  the  axial  and  circumferential 
directions.  The  radial  equilibrium  equation  for  each  ring  and  the  axial  force  equilibrium 
equation  for  the  rotor  are  concurrently  solved  and  expressed  in  a  stiffness  matrix  with  the 
radial  displacements  and  the  uniform  axial  strain.  Once  the  displacements  and  the  axial 
strain  are  obtained,  the  stresses  and  strength  ratios  are  calculated  for  all  the  rings. 


The  stress  distribution  in  each  ring  is  governed  by  the  radial  equilibrium  equation, 
which  is  written  in  cylindrical  coordinates  as  [24] 

d  cr  cr.  cr a  2  n  / 1  \ 

— ^ - Q-  + prco  -  0  (1) 

dr  r 

where  crr  and  <Je  are  the  radial  and  circumferential  stresses,  respectively,  p  denotes  a 
density,  and  co  the  rotational  angular  velocity.  The  stress-strain  relation  in  the  cylindrical 
coordinate  system  as  shown  in  Figure  1  is  written  as 

<T  =  Q(e-aAT)  =  Qs-fiAT  (2) 

where  cr  and  s  are  the  stress  and  strain  vectors,  respectively,  Q  and  a  denotes  a  stiffness 
matrix  and  thermal  expansion  coefficient  vector  in  cylindrical  coordinates.  Eq.  (2)  can  be 


39 


rewritten  in  the  cylindrical  coordinate  system  as 
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In  Eq.  (3),  the  axial  stress  is  zero  in  a  plane  stress  (PSS)  assumption  and  the  axial  strain  is 
zero  in  a  plane  strain  (PSN)  assumption.  More  generally,  the  axial  strain  is  assumed  to  be 
a  constant  in  the  generalized  plane  strain  (GPS)  assumption.  Together  with  these 
assumptions,  the  circumferential  and  radial  strains  are  expressed  by  the  r-direction 


displacement  ur: 


ur  dur 

£e  ~  >  £r  ~~T~’ £& 

r  or 


=  0,  and  ( 


crz  =0  PSS 
sz=  0  PSN 


£=£„  GPS 


(4) 


Substituting  Eqs.  (3)  and  (4)  into  Eq.  (1)  yields  the  radial  displacement  ur\ 


ur  =  - pco 2 (p0r3  +  Cx(pxrK  +  C2(p2r  K  +  cp}£or  +  q>nATr  (5) 


Ci  and  C2  in  Eq.  (5)  are  unknown  constants  to  be  determined  from  the  boundary 
conditions,  and  k  and  <p.  are  defined  in  terms  of  the  material  properties: 
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Substitution  of  Eq.  (5)  into  Eqs.  (3)  and  (4)  yields 

£r  -  -3 poPcpp1  +  /cCx<px rK~'  -  KC2(p2r~*~x  +  <p2£0  +  (pnAT 
or  =  -porcp^r1  +  CxrK~'  +  C2r~K~x  +  (ps£o  +  <pT2AT 


where 


(6) 


(7) 

(8) 


<Pa  =  iP” \Qn  ’  ^=(031+  Q33 ) <P3  +  Q32  and  <pT2  =  (031  +  033 ) (pTX  - pr  (9) 
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In  order  to  derive  the  stiffness  matrix,  the  radial  displacements  at  the  inner  and 
outer  surfaces  are  written  in  a  displacement  vector  using  Eq.  (4), 

u  =  paP'u^  +  £auc  +  A  Tut  +  G0C  ( 1 0) 

where 
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(11) 


In  Eq.  (11),  r,  denotes  the  inner  radius  of  a  ring,  and  ra  the  outer  radius  of  a  ring.  The 
normal  radial  stresses  at  those  two  radii  can  be  also  written  in  a  force  vector  using  Eq.  (8), 

ft  =  P^fo  +  eJi  +  A  Tft  +  1'GC  (12) 


where 
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Eliminating  the  constants  C  in  Eqs.  (10)  and  (12)  derives  the  relation  of  stresses  and 
displacements  at  the  inner  and  outer  surface  of  each  ring.  The  derived  equation  can  be 
expressed  as 


ku  =  fb  +  pa )2fa  +  ejc  +  AT/t 


The  ring  stiffness  matrix  k  is 


k  =  I*G0'G' 
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where 
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Notice  that,  in  defining  the  force  term  ft  of  Eq.  (14),  the  stress  component  was  multiplied 
by  its  radius,  and  the  ring  stiffness  matrix  k  becomes  symmetric.  In  Eq.  (14 ),/®,/f  and  fr 
are  related  to  the  centrifugal  force  and  defined  as 


(  f  \ 
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=  -fa+kUa, 
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(17). 


fr  ~ 


fr, 


=  -ff+kuT 


The  constant  axial  strain  ea  in  Eq.  (14)  can  be  computed  from  force  balance  in 

the  axial  direction.  Since  there  is  no  external  axial  force  applied  on  the  top  and  bottom 
surfaces,  the  summation  of  the  axial  forces  over  all  the  rings  should  be  zero: 


(18) 


Using  the  expressions  of  az  in  Eq.  (3),  and  ur  and  er  in  Eqs.  (5)  and  (7),  Eq.  (18)  is 
expressed  as 


In  Eq.  (19),  j  denotes  a  ring  number  and  q  of  a  y'-th  ring  is  defined  as  follows  (the 
superscript  j  is  omitted  for  clarity): 

q  =  gT0-'G'  (20) 


where 


(21) 


In  Eq.  (19)  and  Eq.  (21),  the  following  symbols  are  defined: 
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Eqs.  (1)-(19)  apply  to  any  arbitrary  ring.  At  each  interface  between  two  adjacent 
rings  (e.g.,  j  and  y+1),  the  continuous  radial  traction  and  continuous  radial  displacement 
should  be  satisfied: 


<7{J*V)  =  aij)  and  u[J+])  =  u[J) 


(23) 


Combining  all  the  ring  stiffness  matrix  in  Eq.  (14)  by  using  the  continuity  equations 
yields  the  global  stiffness  matrix  equation 
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where 


Kuu=tJkU)>Kuc=Kcu  =  fjq(i)T  and  K„  =  fixW -qV?  u?)  (25) 

7= 1  7=1  7=1  V  ’ 

The  global  displacement  vector  U  consists  of  N+ 1  displacements  at  all  the  radii  of  the 
rings.  The  global  force  vectors  Fa  and  FT  in  Eq.  (24)  are  as  follows; 


7=1 
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(26) 


h  i  v  7 

^r  =  I  /r(;>  and  F;  -  £  *W ) 

Notice  that  the  global  stiffness  matrix  is  symmetric,  and  Ft,  has  only  the  boundary 
pressures  given  at  the  inner  and  outer  surfaces  since  the  internal  stresses  were  canceled 
out  due  to  the  stress  continuity  condition.  Fm  is  an  internal  inertial  force  vector  due  to  the 
rotation. 

Now  we  can  solve  the  system  of  N+2  equations  and  obtain  the  radial 
displacements  iP  at  each  interface  and  constant  axial  strain  sa .  The  strains  in  cylindrical 
coordinates  are  then  calculated  using  the  relationship  in  Eq.  (4). 

The  on-axis  strains  e  in  the  material  symmetric  axes  are  obtained  by 
transforming  the  off-axis  strains  e  in  the  cylindrical  coordinates[25]  as  shown  in  Figure 
2: 


e  =  Te 

where  T  is  a  transformation  matrix  defined  in  terms  of  a  ply  angle  (f> ; 

cosV  sin2^  0  cos^sin^ 

sin2^  cos 1  <j>  0  -cos^sin^ 

0  0  1  0 

-2  cos  ^  sin  ^  2  cos  (f> sin  (j>  0  cos2  ^-sin2  (j> 

The  on-axis  stress  a  is  then  calculated  using  the  stress-strain  relationship: 

<x  =  Q(s  -a AT) 


(27) 


(28) 


(29) 


The  on-axis  stiffness  matrix  Q  and  the  coefficients  of  thermal  expansion  (CTE)  a  in  Eq. 
(29)  are  given  by  the  material  properties[25]: 
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The  definitions  of  Moduli,  Poisson’s  Ratios,  and  CTE  in  Eq.  (30)  are  explained  in  Table  1. 

Once  the  on-axis  stresses  in  the  y-th  ring  are  calculated,  a  failure  of  the  material 
can  be  accessed  by  using  a  failure  criterion.  In  this  study,  3-dimensional  Tsai-Wu 
quadratic  failure  criterion[19]  is  used: 

aTFa  +  FaR-R2  =  0  (31) 


Solving  Eq.  (31)  yields  a  strength  ratio  R  indicating  material  failure  with  a  value  greater 
than  1 .  When  the  strength  ratio  R  <  1 ,  the  applied  stresses  can  increase  by  a  factor  of  its 

inverse  before  material  failure  occurs.  The  strength  parameters  F  and  F  are  defined  in 
terms  of  the  material  strengths: 


where  X  and  X '  are  the  tensile  and  compressive  strengths  in  the  fiber  direction,  Y  and  Y' 
those  in  the  matrix  direction,  and  S  the  shear  strength. 
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Figure  1.  Model  of  multi-ringed  composite  flywheel  rotor  with  the  angle  ply  in  each 

rings.  The  radial  displacement  and  stress  of  j~th  element  components,  and  the  off-axis 
coordinate  (0,  z,  r)  and  the  on-axis  coordinate  (  0,  z,  r)  systems  are  defined. 
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2.  OPTIMIZATION  COMPUTER  CODE  -  FLYOPT 


A  proper  combination  of  composite  materials  for  the  rotor  could  reduce  the  stress 
distribution  in  the  rotor  and  eventually  increase  TSE.  The  objective  of  this  study  is  to 
maximize  TSE  of  a  hybrid  flywheel  rotor  without  any  material  failure.  Strength  ratio  R  is 
used  for  failure  analysis,  where  a  value  of  R  greater  than  1  indicates  material  failure  [8]. 
To  ensure  the  rotor  against  any  material  failure,  in  this  approach,  strength  ratio  R  and  the 
damage  indices  Z^and  /^calculated  at  any  point  of  the  rotor  must  have  values  less  than 

1 .  The  inner  radius  n  and  height  h  of  the  rotor  are  given  and  the  thickness  tm  of  each  rim 
(m=l,2,3,...,A/)  is  used  as  the  design  variables.  The  optimization  is  then  formulated  as 
follows: 


Find  tm  (m  =  1,2, 3, 

Maximize  E  -—Ico2  (50) 

2 

Subject  to  Max  R,Da,De<  1 

where  I  is  a  mass  moment  of  inertia  of  the  hybrid  rotor,  and  co  is  a  rotating  speed.  The 
total  stored  energy  (TSE)  E  can  be  thus  written  in  terms  of  the  inner  radius  rm  and  outer 
radius  rm+i  of  each  m-th  rim: 


E  = 


(51) 


where  the  subscript  m  indicates  each  rim  and  the  other  variables  include  the  length  of  the 
rotor  h,  and  the  density  of  the  m-th  rim  pm . 

In  order  to  solve  the  nonlinear  optimization  problem  formulated  in  Eq.  (50),  the 
modified  method  of  feasible  directions  for  constrained  minimization  is  used.  The  main 
task  in  the  optimization  method  is  to  find  a  usable-feasible  search  direction  and  the  one¬ 
dimensional  search.  In  the  every  iteration  of  optimization,  stress  and  strength  analysis  are 
performed  and  the  maximum  strength  ratio  in  each  ring  is  evaluated  as  described  below. 
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An  optimal  design  of  a  hybrid  composite  rotor,  by  maximizing  TSE  without  any 
material  failure  as  formulated  in  Eq.  (50),  is  performed  using  a  combination  of  composite 
materials  with  different  fiber  types:  (A)  Glass/Epoxy,  (B)  T300/2500,  (C)  T800H/2500. 
Optimizations  are  performed  combining  up  to  three  different  composite  rims  and  can  thus 
be  grouped  into  3  categories.  In  the  first  category,  only  one  material  is  used,  and  in  the 
second  category,  two  materials,  and  so  on.  In  each  category  only  several  cases  of  material 
selections  and  sequences  are  considered  because  it  is  efficient  way  of  increasing  TSE  to 
arrange  softer  materials  in  the  inner  side  and  stiffer  materials  in  the  outer  side  of  the  rotor. 
The  thickness  of  each  composite  rim  is  considered  as  a  design  variable  and  equally 
divided  into  element  rings,  resulting  in  a  total  of  30  element  rings.  The  initial  thickness 
for  each  composite  rim  is  20  mm,  and  the  inner  radius  of  the  rotor  is  50mm.  The  value  for 
the  rotating  speed  is  taken  as  60,000  RPM. 

All  the  calculated  optimum  TSE  are  shown  in  Figure  2  where  each  bar  indicates 
the  location  of  each  rim  and  thus  the  length  of  the  bar  means  the  thickness  of  each  rim. 
As  noted  in  the  stress  distribution,  the  effects  of  the  curing  temperature  on  the  optimal 
design  are  significant.  For  instance,  the  results  of  the  optimal  designs  using  the  material 
sequence  A-B-C  show  that  TSE  reduces  by  about  30%  with  AT  =  -110°C  and  36%  with 
AT  =  -125  °C.  In  other  words,  neglecting  the  residual  stresses  may  leads  to  causing  the 
higher  stresses  than  expected  in  the  process  of  operating  the  flywheel  rotor.  The 
distributions  of  the  strength  ratios  for  A T=  0°C,  AT  =  -110°C  and  AT  =  -125  °C  are 
shown  in  Figure  6.  It  is  noted  that  the  maximum  TSE  is  obtained  when  the  strength  ratio 
R=l  is  reached  in  every  composite  rim.  Thus,  strength  is  not  wasted  in  each  rim  which 
otherwise  experience  proportionately  smaller  stress  levels.  This  means  that  maximizing 
TSE  is  equivalent  to  incorporating  the  highest  possible  stresses  that  satisfy  the 
equilibrium  and  compatibility  conditions.  At  the  optimum  state,  each  thickness  of  the  rim 
is  thus  determined  so  that  the  strength  ratio  R=1  is  reached  in  every  composite  rim.  It  can 
be  obviously  noted  that  the  number  and  sequence  of  materials  are  very  important  in 
maximizing  TSE  or  minimizing  the  stresses  inside  the  materials.  The  stress  distributions 
are  effected  by  the  sequence  of  stiffness,  mostly  hoop  directional  Young’s  modulus. 
Softer  materials  inside  of  the  rotor,  which  are  easier  to  expand  in  the  radial  direction,  and 
stiffer  materials  outside  of  the  rotor,  which  experience  less  radial  displacement, 
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eventually  reduce  the  radial  stresses. 

A  computer  program  FLYOPT  is  developed  using  Microsoft  Excel  user-interface. 
The  input  data  in  blue  and  the  output  data  in  red  for  the  user-interface  for  the 
optimization  are  respectively  shown  in  Figures  3.  The  radial  displacement  and  the 
stresses  are  automatically  updated  according  to  the  results  calculated  from  FlyOpt.  The 
design  variables  in  the  optimization  module  can  be  linked  to  any  variables  in  the  stress 
analysis  module. 
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Figure  2.  Comparison  of  optimal  thickness  of  composite  rims  and  the  corresponding  TSE 
(Wh)  for  various  selections  of  composite  material  (iy=60,000RPM  ). 
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Figure  4.  The  distribution  of  the  radial  displacement,  the  stresses,  and  the  strength 
ratios;  The  data  are  automatically  updated  according  to  the  results  from  FLYOPT. 
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3  CONCLUSIONS 


The  goal  of  the  present  work  is  to  develop  an  optimization  tool  for  a  design  of  the 
flywheel  rotor  with  a  consideration  of  the  process-induced  residual  strains  accumulated 
during  the  winding  and  curing  process.  A  stress  analysis  of  a  helically  wound  hybrid 
multi-ring  rotor  is  developed  in  which  an  equivalent  temperature  change  is  used  to 
consider  the  residual  stresses  and  strains.  The  validity  of  the  stress  analysis  has  been 
proved  using  a  three-dimensional  finite  element  program.  The  effects  of  the  material 
hybridization  and  the  helical  winding  angles  on  the  thermal  stresses  are  also  investigated. 
The  residual  stresses  have  been  experimentally  measured  using  a  split-ring  test.  As  a 
product  of  this  work,  a  computer  program  for  the  optimum  design  of  the  hybrid 
composite  flywheel  rotor  are  developed  which  covers  the  various  size  and  material 
properties  of  the  rotor. 
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